Mode-Field Radius of Photonic Crystal Fibers Expressed by the V— parameter 
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We numerically calculate the equivalent mode-field radius of the fundamental mode in a photonic 
crystal fiber (PCF) and show that this is a function of the V-parameter only and not the relative 
hole size. This dependency is similar to what is found for graded-index standard fibers and we 
furthermore show that the relation for the PCF can be excellently approximated with the same 
general mathematical expression. This is to our knowledge the first semi-analytical description of 
the mode-field radius of a PCF. 
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Theoretical descriptions of Photonic Crystal Fibers 
(PCFs) have traditionally been based on numerical meth- 
ods such as the plane- wave expansion method, 0,0 meth- 
ods employing localized functions, 0, Q or the multipole 
method. 5, 6] A numerical approach is generally required 
due to the complex dielectric cross section of the PCF 
which makes analytical approaches very difficult and re- 
sults in the fact that no close-form analytical descrip- 
tions of either propagation constants or mode fields are 
available. In this Letter we investigate the relation be- 
tween the equivalent mode-field radius of the fundamen- 
tal mode and a recently proposed formulation of the V- 
parameter for a PCF. y| We show that the mode-field 
radius is a function of the V-parameter only and provide 
an empirical expression describing this relation having 
the same mathematical form as known from the descrip- 
tion of graded-index standard fibers. @ 

The PCF analyzed in this work is an all-silica fiber 
having a triangular arrangement of circular voids with di- 
ameter d running along the full length of the fiber. The 
voids are placed symmetrically around a central defect 
which acts as the fiber core consisting of a solid silica re- 
gion i.e. an omitted air hole. The air-hole matrix, which 
has a lattice constant, A, functions as a cladding region 
and the fiber structure is invariant in the longitudinal di- 
rection. Such a fiber was first proposed in,9j] and it was 
shown that it posses u niqu e properties such as endlessly 
single-mode operation. 10] The endlessly single-mode op- 
eration is a consequence of the fact that the number of 
guided modes is finite regardless of wavelength and that 
the upper limit for this number decreases with the air- 
filling fraction of the structure i.e. with the value of d/A. 
Sufficiently small air holes will cause the number of al- 
lowed modes to equal two, namely the two deg enerate 
polarization states of the fundamental mode.yjj] 

In the case of graded-index standard fibers, of which 
the step-index fiber (SIF) can be considered a special 
case, the V-parameter plays a central role in the descrip- 
tion of the number of guided modes, the cut-off criterion 



[11] as well as the mode- field radius. In attempt to ob- 
tain an expression for a V-parameter adequate for PCFs, 
approximations based on an equivalent SIF has gener- 
ally been employed [10,12]. Although this approximation 
contains the overall correct physics it fails to describe 
the cut-off properties and has difficulty when it comes to 
determining an appropriate equivalent core radius. Re- 
cently, we proposed a definition of the V-parameter for 
a PCF,Q Vpcf, rejecting the SIF approximation [See 
Ref. for a more detailed discussion on Eq. (1)]: 



Vpcf(A) = ^A^? (A)-n=(A) 



(1) 



In Eq. (1), A is the free space wavelength and n co (A) 
and n c i(A) are the effective indices of the fundamen- 
tal mode and the first eigenmode in the perfect infinite 
cladding structure, respectively. The latter is often de- 
noted the fundamental space-filling mode and can be in- 
terpreted as the refractive index of the cladding material 
in the absence of the core. ^3 Using this definition it can 
be showen that the condition for the higher-order mode 
cut-off can be formulated as Vpcf = tt-Q Both effective 
indices in Eq. (1) are strongly wavelength dependent and 
can not be approximated by constants as in the case of 
standard fibers. It is due to the strong dispersion of the 
effective indices that the PCF can be attributed many 
of its unique properties. In Fig. 1, Vpcf calculated us- 
ing the plane-wave expansion method|l3( is plotted as 
function of the normalized frequency, A/A, for values of 
d/A ranging from 0.30 to 0.70 in steps of 0.05. The hor- 
izontal dashed line in the plot indicates Vpcf = tt and 
thereby the single-mode boundary. Since the variation 
of refractive index of silica, n s ii; ca , over the transparent 
wavelength region is in the order of a few %, the effect of 
material dispersion will only result in a small shift in the 
effective indices n co (X) and n c i(A). Furthermore, such a 
perturbation will shift both further reducing the effect on 
VPCF. In the calculations we therefore used a fixed value 
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FIG. 1: Calculated Vpcf (from Eq. 1) as function of the 
normalized frequency, A/A, for d/A equal to 0.30, 0.35, 0.40, 
0.45, 0.50, 0.55, 0.60, 0.65, and 0.70 from bottom and up. 
The dash line indicates the single-mode boundary, Vpcf = 7r. 



of n s iii C a = 1.444, thereby preserving the scale invariancc 
of Maxwell's equations. 

Vpcf approaches a constant value, dependent on d/A 
for increasing A/A|10| and since the number of modes 
generally increases with the ^-parameter, the asymp- 
totic behavior of Vpcf is consistent with the end- 
lessly single-mode property. In the framework of stan- 
dard fibers, the weak wavelength dependency of the V- 
parameter will lead to a mode-field radius which is also 
only weakly dependent on the wavelength. 

For a graded-index standard fiber, the index profile 
can be described by a power law in which case the shape 
of the index profile is governed by an exponent, g. For 
g = 1 the index profile is triangular, for g = 2 the pro- 
file is parabolic, and in the limit of large values of g the 
profile approaches that of a SIF. For this type of fiber, 
the mode field is generally very close to a Gaussian dis- 
tribution provided that the field does not penetrate too 
deeply into the cladding region. The agreement between 
the actual mode and a Gaussian distribution has a weak 
dependency on g and is a perfect match in the case of 
a parabolic index profile. The mode-field radius, w, is a 
function of the V-parameter and can be fitted using the 
expression: j| 
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A, £>, and C, are fitting parameters which are depen- 
dent on the index profile while a denotes the core radius. 
In the case of a SIF the first term in Eq. (2) is constant 
and the values of A, B, and C are 0.65, 1.619, and 2.879, 
respectively. The accuracy of the fits provided by Eq. 
(2) is better than 2% in the range 1.5 < V < 7 and the 
relation is therefore extremely useful when designing and 
working with SIFs or graded index fibers in general. 



In order to investigate a similar relation for the PCF, 
we first introduce an equivalent mode-field radius of the 
fundamental mode, wpcf, as the mode-field radius of 
the Gaussian distribution having the same effective area, 
A e s, as the fundamental mode itself, yielding the simple 
relation A c g — ttw 2 . The effective area is calculated as 



A, 



dr i I(r 
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(3) 



In Eq. (3), I{r±) is the transverse intensity distribu- 
tion of the mode. Although the intensity distribution of 
the fundamental mode in a PCF is not rotational sym- 
metric but rather has the 6-fold symmetry of the tri- 
angular cladding structure, a Gaussian approximation is 
in fact very good and has previously been employed in 
the description of various PCF properties. [H El El By 
numerical calculation of A e g as function of the normal- 
ized wavelength, A/A, the normalized mode-field radius, 
^pcf/A can be plotted as function of Vpcf as shown 
in Fig. 2. The normalization with A is chosen since A 
is the natural length scale of the problem in the same 
sense as the core radius, a, is for the graded-index fiber. 
The plot in Fig. 2 represents data for d/A — 0.30 and 
d/A= 0.70 indicated by open squares and circles, respec- 
tively. The data points for these two calculations overlap 
in the entire range where the value of Vpcf overlap and 
data for calculations having d/A — 0.35, 0.40, 0.45, 0.50, 
0.55, 0.60, and 0.65 (not shown) fall on the same curve 
indicating that wpcf /A is in fact a function of the Vpcf 
only. Also included in Fig. 2 is the corresponding curve 
for the SIF (dashed line) showing wsif/o as function of 
Vsif calculated from Eq. (2). 

The functional dependency of w pcf /A on Vpcf is seen 
to follow the same overall behavior as wsif/o as function 
of Vsif- For increasing values of Vpcf, wpcf/A must 
saturate at a constant value corresponding to the mini- 
mum obtainable mode size whereas it grows dramatically 
for sufficiently small values of Vpcf indicating that the 
mode is weakly confined to the core and begins to pene- 
trate into the cladding region. This is in good agreement 
with the fact that small values of Vpcf are realized ei- 
ther for small values of d/ A or if the optical wavelength 
is similar to the structural dimensions, A ~ A. 

Due to the similar functional dependency of the mode- 
field radii on the respective ^-parameters it is reasonable 
to use the same fitting function for the PCF as used for 
graded-index standard fibers. However, if assuming the 
simplest case of a SIF it is not possible to obtain a good fit 
and the more general expression, for which g can assume 
any positive value, was therefore employed. This results 
in the parameters A = 0.7078, B = 0.2997, and C = 
0.0037, respectively in the case for g = 8. The fit based 
on these parameters are plotted in Fig. 2 (solid line) 
and is seen to be very good for small as well as for large 
values of Vpcf- The maximal deviation between values 
predicted by the fit and the data points is determined to 
be less than 1%. 
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FIG. 2: Normalized mode- field radius for the PCF, wpcf/A, 
as function of Vpcf for d/A — 0.30, (squares) and d/A — 0.70 
(circles) along with a numerical fit through these points (full 
line) . The corresponding expression for the normalized mode- 
field radius, wsw/a for the SIF as function of Vsif is also 
shown (dashed line). The two arrows indicate the higher- 
order mode cut off of the SIF and PCF at Vsif = 2.405 and 
Vpcf = f, respectively. 



An important difference between the PCF and the SIF 



can be learned from the curves in Fig. 2 on which the 
higher-order mode cut-off are indicated at Vpcf = T an d 
Vsif = 2.405, respectively. Whereas the cut-off for the 
PCF is located at a point where the curve is relatively flat 
the SIF cut off is located where the curve is much steeper. 
Also, Vsif depends on the wavelength as A -1 further in- 
creasing the wavelength dependency of the mode size in 
the single-mode region of the SIF. For the PCF, not only 
the mode-size dependence on Vpcf close to cut off is week 
but also the dependence of Vpcf itself on A provided that 
A is a few times larger than A or more (See Fig. 1). This 
results in a mode- field size which is essentially wavelength 
independent for the PCF compared to the SIF. 

In conclusion, we have shown that the recently pro- 
posed expression for the V-parameter of a PCF uniquely 
determines the normalized mode-field radius. Further- 
more, we have presented a simple analytical expression 
describing this relation based on the same general expres- 
sion used for graded-index standard fibers and thereby 
further closed the gap between standard fibers and PCF 
in terms of available tools describing their properties. 



M.D. Nielsen acknowledges financial support by the 
Danish Academy of Technical Sciences. M.D. Nielsen's 
e-mail address is mdn@crystal-fibre.com. 



[1] J. Broeng, S. E. Barkou, A. Bjarklev, J. C. Knight, T. A. 
Birks, and P. St. J. Russel, Opt. Commun. 156, 240 
(1998). 

[2] A. Ferrando, E. Silvestre, J. J. Miret, and P. Andres, 

Opt. Lett. 24, 276 (1999). 
[3] D. Mogilevtsev, T. A. Birks, and P. St. J. Russel, Opt. 

Lett. 23, 1662 (1998). 
[4] T. M. Monro, D. J. Richardson, N. G. R. Broderick, and 

P. J. Bennett, J. Lightwave Technol. 17, 1093 (1999). 
[5] T. P. White, B. T. Kuhlmey, R. C. McPhedran, D. 

Maystre, G. Rcnversez, C. M. de Sterke, and L. C. Bot- 

ten, J. Opt. Soc. Am. B 19, 2322 (2002). 
[6] B. T. Kuhlmey, T. P. White, G. Renversez, D. Maystre, 

L. C. Botten, C. M. de Sterke, and R. C. McPhedran, J. 

Opt. Soc. Am. B 19, 2331 (2002). 
[7] N. A. Mortensen, J. R. Folkenberg, M. D. Nielsen, and 

K. P. Hansen, Opt. Lett. 28, 1879 (2003). 
[8] D. Marcuse, J. Opt. Soc. Am. 68, 103 (1978). 



[9] J. C. Knight, T. A. Birks, P. St. J. Russel, and D. M. 
Atkin, Opt. Lett. 21, 1547 (1996). 
[10] T. A. Birks, J. C. Knight, and P. St. J. Russel, Opt. Lett. 
22, 961 (1997). 

[11] W. A. Gambling, D. N. Payne, and H. Matsumyra, Elec- 
tron. Lett. 13, 139 (1977). 

[12] J. C. Knight, T. A. Birks, and P. St. J. Russel, J. Opt. 
Soc. Am. A 15, 748 (1998). 

[13] S. G. Johnson and J. D. Joannopoulos, Opt. Express 8, 
173 (2001). 

[14] N. A. Mortensen and J. R. Folkenberg, Opt. Express 10, 
475 (2002). 

[15] N. A. Mortensen, Opt. Express 10, 342 (2002). 

[16] N. A. Mortensen, J. R. Folkenberg, P. M. W. Skovgaard, 

and J. Broeng, IEEE Photonic. Tech. Lett. 14, 1094 

(2002). 



